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Abstract. We develop a local moment approach to static properties of the symmetric 

Anderson model in the presence of a magnetic field, focussing in particular on the strong 
coupling Kondo regime. The approach is innately simple and physically transparent; but is 
found to give good agreement, for essentially all field strengths, with exact results for the 
Wilson ratio, impurity magnetization, spin susceptibility and related properties. 



JZZ \ Theoretical study of Anderson impurity models (AIM) M has long formed a cornerstone 

of research in strongly correlated electron systems |2J . Its pursuit has recently acquired added 

£SJ , impetus both from the experimental discovery that direct, tunable mesoscopic realisations of 

C^- ■ AIMs are provided by nanostructure devices such as quantum dots []3[ , or surface atoms probed 
by scanning tunneling microscopy [|4|]; as well as from the advent of dynamical mean- field 
theory for correlated lattice-based fermions ||,[|, within which models such as the Hubbard 
or periodic Anderson lattices map onto an effective, self-consistent AIM. There is however a 
disparity between progress in theoretical descriptions of static (thermodynamic and related) 
properties of AIMs, and dynamical properties such as single-particle excitations. The former 
are generally well understood ||, at least for the conventional (metallic) AIM uj], using a 
variety of powerful techniques that often yield exact results but unfortunately do not in 
general permit ready access to dynamics. Extant theories for the latter by contrast, which 
in principle should also recover static properties as a limiting case and are approximate by 
necessity, suffer from well known qualitative limitations gj. And while numerical approaches, 
notably Wilson's numerical renormalization group (NRG) [M, circumvent most beautifully this 

j^ ■ impasse and provide benchmark numerical results (see e.g. [0), the need for new, approximate 
H \ theoretical methods is not thereby vitiated. 

We have recently initiated development of one such — the local moment approach (LMA) 
P, f| - a non-perturbative many-body method in which the notion of local moments [fjj 
is introduced explicitly and self-consistently from the outset. The underlying approach is 
simple and physically transparent, but not without virtue: for example it handles single- 
particle dynamics on all energy scales, while recovering correctly Fermi liquid behaviour at 
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low energies; capturing in particular the spin-fluctuation physics charactristic of the strong 
coupling Kondo regime, where it leads to good quantitative agreement with NRG results 
for e.g. the universal single-particle scaling spectrum (see fig. 10 of PQ] ), The LMA is not 
moreover specific to the Fermi liquid behaviour ubiquitous in the conventional AIM. It can 
also handle non-Fermi liquid phases such as arise in the soft-gap AIM, to which the approach 
has also recently been applied M and compared successfully to NRG results in mn . 

Static properties resultant from the LMA have not however been considered hitherto. That 
is the purpose of the present paper: to consider the LMA for the conventional symmetric 
AIM, generalised to include the effect of an applied magnetic field, in order to determine 
static properties and their field dependence, notably the excess impurity magnetization (and 
hence linear conductance of a quantum dot (see below)), the spin susceptibility and Wilson 
ratio. Our primary emphasis is naturally on the strong coupling Kondo regime, where we 
make comparison to and and find striking agreement with exact results provided by the Bethe 
ansatz pi| . This provides an acid test for any approximate theory, and since the LMA is 
confined neither to statics nor to integrable models, its success in describing quantum impurity 
models for which exact results are available provides important evidence for the strength of 
this new approach in generic situations where exact results will not typically be possible. 

The AIM Hamiltonian is given in standard notation by 

H = ^ € k fl k<J + ^ Uia + ^i-o) flier + ^ V ik Ul a C kcr + Il.C.J . (1) 

k.a a k.a 

The first term refers to the non-interacting host; and the second to the impurity with on-site 
interaction U. ti„ — e* — ah includes the (local) Zeeman coupling to the external field H, 
with h = ^gfisH and a = +/— for f / |-spin electrons. For the symmetric AIM considered, 
ei = — ¥; and by particle-hole (p-h) symmetry the impurity charge u, = y~] < "hia >= 1 V /i. 
The final term in eq. (fil) is the host-impurity coupling. The total impurity Green function 
G(ui; h) = \ 2 G a (uj, h), with G a (uj; h) — G^ (w; h) — i sgn(w)7rD cr (a;; h) given by 

a 

G CT (w; h) = \lu+ - A(w) + ah - E CT (w; h)] (2) 

where uj + = ui + i0 + sgn(w). Here A(u>) = Ar(cj) — i sgn(tu)Ai(u>) is the host-impurity hy- 
bridization, with Ai(u>) — it ^2 k V^ k 5(u! — efe); in practice we consider explicitly the wide-band 
AIM with Aj(u) =A Vw (and A R (w) = 0). t a {w;h) = T%(w;h) - i sgn(cj)tl(cj; h) de- 
notes the self-energy (excluding the trivial Hartree contribution); by p-h symmetry E,j(cl;; h) = 
— E_ CT (— uj\ h) and D a (uj;h) = D^ a (—uj;h), such that D CT (0;/i) = D(0;h) at the Fermi level 
lu = 0. The h-dependent (but a- independent) quasiparticle weight is defined by z(h) = 

[l-(flE*teW^)^~\ 

While H a (u>;h) in its entirety naturally determines spectral dynamics, its u) — limit 
enables the /i-dependence of a range of important static properties to be deduced. The single- 
particle spectrum at the Fermi level is given by ttAqD(0; h) = [l + ([h — aY^(0; /i)]/Ao) ] ; 
it determines the (T = 0) linear conductance of a quantum dot modelled by an AIM, via JE2J 
G c (h) = (e 2 /irh)irAoD(0;h) (for equal dot-lead couplings). Likewise the (excess) impurity 
magnetization M[(h) follows from the Friedel sum rule as 0| 

Mi (h) - ^tan- 1 ([h- at*(0;h)] /Ao) (3) 
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(for either spin a), with the corresponding susceptibility Xi{h) = dMJdH given by 



X i(h) = i( 5MB ) 2 £(0; h) [1 - <r(dl%(p; h)/dh)] . 



(4) 



And the Wilson ratio, Rw(h) = axi(h) / ji(h) with ji(h) the linear specific heat coefficient 

(a = [2Trk B ] 2 /[3(gfi B ) 2 }), is given by 



R w (h) = z(h) [1 - <r(9E*(0; h)/dh)] . 



(5) 




Within the LMA |,|], E*(w;/i) = -§L/|/x(/i)| + £ CT (w;/i) 
is separated into a purely static Fock contribution (with local 
moment |/j(/i)|) that alone would survive at the crude mean- 
field (MF) level of unrestricted Hartree-Fock; plus a dynamical 
contribution H a (uj)h). The latter includes in particular a non- 
perturbative class of diagrams (fig. Ill) that embody dynamical 
coupling of single-particle excitations to low-energy transverse 
spin fluctuations, and hence capture the spin-flip scattering es- 
sential to describe the strong coupling (SC) or Kondo regime for 
U = U/ttAq ^> 1. These are expressed in terms of MF propaga- 
tors (solid lines in fig. pi), 



Fig. 1 - LMA E<r(w), see 
text. Wavy lines denote U. 



g a {uo-h)= uj+ - A(w) + a(~U\n(h)\ + h) 
and for ui = in particular £j(0; h) ( = S?(0; h)) is given by 



(G) 



£f(0;/i) = C/ 2 



divi . 



Imn+-(wi) e(u>i)Reg7{ojf,h) +6(-u>i)Regf(u)i;h) 



(7) 



The transverse spin polarization propagator II + (u>) (shown hatched in fig. |l|) is given at the 
simplest level by an RPA-like particle- hole ladder sum in the transverse spin channel; Go- (cu) 
denote the one-sided Hubert transforms of Q a {u)\ h). 

LMA : h — 0. We first review briefly the LMA for h = 0, an understanding of which 
underpins the general case ft. ^ 0; full details are given in H. Here the key idea is symmetry 
restoration, reflected in S?(0; 0) — S?(0; 0); i.e. (using p-h symmetry) 



E?(0;0) = --ff| M (0)|+E?(0;0) = 0. 



(8) 



Imposition of eq. (|3|) as a self-consistency condition, achieved in practice for given U by 
varying the local moment |^(0)| from its MF value, preserves the [/-independent pinning of 
the Fermi level spectrum (7rAo-D(0; 0) = 1); and in turn leads to Fermi liquid behaviour at 
low-energies. Most importantly, it introduces naturally a low-energy spin-flip scale 10^ that 
sets the timescale for restoration of the broken symmetry endemic at MF level, and is manifest 
in a strong resonance centred onw = cj^ in Imn -1 (lu). This is the Kondo scale. Its form in 
SC is readily deduced via eq. (^), which for U ^> 1 has the asymptotic form || 
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S R (0 ;0) = ^ln 

■ TT 



u_ 

IT 



(9) 



Combining with eq. (|S|) (using |/x| — » 1 in SC) gives the SC Kondo scale for h — 0, viz 
w^ ~ U exp ( — 7rC//8Ao) where the exponent is exact (the prefactor, in effect a uv-cutoff, is 
of course approximate). The leading low-w behaviour of Ej* (w; 0) is also readily obtained f| 
in SC: (9S^(w;0)/9w) w ~ — 4Ao/7rw° 1 , producing a quasiparticle weight z(0) ~ ttuj^/AAo 
that is naturally proportional to the Kondo scale. 

LMA: h > 0. For /i = 0, two degenerate MF states arise [g|. The degeneracy is naturally 
removed for h ^ 0, where one or other state is picked out according to sgn(/i). The appro- 
priate MF propagators then determine the many-body LMA self-energies Y, a: which satisfy 
Y, a (uj;h) = £_ ff ((i>; — h) as well as p-h symmetry (and hence e.g. M;(/i) = —Mi(—h)). For 
convenience in what follows we consider explicitly h > 0. Our primary aim is to consider the 
SC Kondo limit of the AIM. To this end we thus consider finite h/A o z(0) in the SC limit 
where z(0) — ► 0. 

Since G a — [G^ 1 — E ff ] , and the LMA S^ = E \Qa\ is a functional of the MF propaga- 
tors given by eq. (g), it follows generally that the h-dependence is embodied fully in x(h) = 
f \fj,(h)\ + h (with |/i(0)| such that eq. (g) is satisfied for h = 0). The LMA £ CT (u>; h) for h > 
is thus formally equivalent to that for /i = 0, but with x(0) — ► z(/i) = #(0) + [^5|/a(/i)| + h] 
(where 6\n(h)\ — \fi(h)\ — |/x(0)| is calculated in practice at MF level). In consequence, S^(0; h) 
and (9E^-(w; h)/du>) in SC have idential forms as for h = 0, but with w^ = w m (ft = 0) 
replaced by w m (/i) = w m (x(/i)) where U) m (h) is again the spin-flip scale characteristic of 

Imn + -(w); i.e. Ef(0;fc) ~ ^f 2 - In [U/w m (hj\ , and (dE*(u;; fc)/0w) o 4A /mj m (/i) such 

that z(ft,)/z(0) = W m (/i) / 'w^. Hence, using eq. (||) for h — 0, 



TT 



*(fc) 



z(0) 



(10a) 



For finite h/Aoz(0), the — -^(5|/j,(/i) | contribution to S?(0; h) is readily shown to be exponen- 
tially small in SC and is thus neglected in eq. (10a). By the same token uj m (h) = ui m (x(h)) 
satisfies dw m /d/i = du) m /dx; with duj m /dx — 2 in SC known from the LMA for h — S. 



Hence . 



>w 



2/i, 



*(0) 



7T 
1 + -- 



2 A o z(0) 



feB^b 



(10b) 



Here we have also introduced the scale Tq defined by fee To = 7rAo2(0)/4, in terms of which 
the exact linear susceptibility in the SC Kondo limit is ||,|l|] Xi(0) = (5Mb) 2 /4£;bTo. 

Eqs. (10) are the basic LMA equations appropriate to the SC Kondo limit of the AIM. 
The Wilson ratio follows directly from eqs. (|S|) and (10). It is given by Rw(h) = 2 for all 
finite h = h/k^To and feeT) — > 0. This is the exact result for the Kondo limit, including its 
field (in)dependence J13[ |. It is not however recovered by conventional approaches such as the 
slave boson mean-field approximation or -^-expansions to leading order (with N the impurity 
degeneracy), which give M i?\y(0) = 1 corresponding to non- interacting quasiparticles. 

The magnetization in SC likewise follows from eqs. (10) with (0) (where the 'bare' h may 
be neglected since h — hk B T — ► 0); namely (for h > 0) 
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Mi (h) = ^ tan" 

7T 



- In | I 

7T 



fe^O 



(11) 



This is exact asymptotically for both weak and strong fields h = h/ksTQ. For h <C 1 it yields 
Mj(h) ~ Xi(0)J? with Xi(0) = (gMB) 2 /4fc B To, thus capturing the exact Xi(0)- For h >■ 1, 
cq. (|ll|) gives 



M; (ft) 



2 



1 



2 In ^fr 



0([lnft] 



(12) 



It thus recovers the slow logarithmic approach to saturation, although only the leading such 
correction is obtained, the exact M\(h) being given by |ll| 2Mjg^ = 1 — [21n(y / 7re/i)] _1 + 
0(lnln/i/[ln/i] 2 ). Fig. || shows the full /i-dependence of M\(h) compared to the exact Bethe 
ansatz result for the Kondo model p3|. The agreement is very good throughout; the deviation 
is at most ~ 10%, occuring for fields h ~ 0(1). Hence the field-dependence of the linear 
conductance, given (using eq. (ph) by G c (h) = (e 2 /irTi) cos 2 [it Mi{h)/g[iB\ , is in turn well 
captured by the LMA, recovering the exact limits (irU / e 2 )G c (h) ~ 1 — [irh/2] 2 for h <C 1 
and ~ [-ln(ft,)]~ 2 for h ^§> 1. Similar comments apply to the susceptibility Xi{h) which is 
also shown in fig. 0, and is likewise exact asymptotically in the weak and strong field limits. 
Finally, since the LMA recovers correctly Rw(h) — 2 for all h, the /i-dependence of the specific 
heat coefficient 7i(/i) oc Xi(/i) is precisely that of XiC 1 )- 
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Fig. 2 - Mi(h)/g(iB and Xi(ft) = [4>kBT /(gHjsf]Xi(h) vs h = h/k B T in SC Kondo limit. Solid line: 
LMA. Dashed line: exact Bethe ansatz [|Ll| for Kondo model. 
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Having discussed the SC Kondo limit, we now turn briefly to LMA results for fmite-t/, 
focusing on the susceptibility Xi{h) = Xi(h;U) (eq. (||)). Fig. || shows the h — h/k^T^ 
dependence of Xi{h;tJ) — [AkBT /{g^) 2 ]xi{h]U) for U — 10, 4 and 0; as above k^T a = 
7rAo2;(0)/4, with z(0) the [/-dependent, zero-field quasiparticle weight. For the trivial non- 
interacting limit (where z(0) = 1), Xi(/i;0) = i[l + (irh/i) 2 }^ 1 and \h = -^-. Note first 

that for h = 0, Xi(0; U) = iR w (0) follows generally from eqs. @), (|) using ttA o D{0, 0) = 1. 
The LMA recovers correctly i?w(0) = 2 in the Kondo limit, but does not yield the exact 
exponential approach (in U) to it, giving instead algebraic behaviour (i?w(0) ~ 2 — 0(l/t/)); 
Xi(0; U) is thus underestimated somewhat by the LMA. Nonetheless for U — 10, and over 
essentially the entire /i-range shown in fig. ||, Xi 1S very close to its universal SC Kondo limit 
(fig. ||, solid line), exhibiting in particular the high-field behaviour Xi ~ [Mn 2 (/i)] _1 for ht 10 2 . 
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Fig. 3 - Xi(h) = [4fc B ro/(g/i B ) 2 ]Xi(' 1 ) v s h/k B T for: the Kondo limit (solid line), U = 10 (point-dash), 
4 (dashed), and (dotted). For U = 4, h = U/2 is marked by a vertical arrow. 



For sufficiently large fields at any finite U however, Xi(h',U) should exhibit a crossover 
to behaviour charactristic of the free-orbital regime. That the LMA correctly captures such 
behaviour is illustrated in fig. [| for U = 4, where h = §■ is marked by an arrow. X;(/i; U) is 
seen to cross over from Kondo-like behaviour and for ft> | becomes proportional (parallel in 
fig. H) to the non-interacting Xi(h;0) (shown dotted in fig. g). This is indeed the free-orbital 
regime: here the susceptibility itself, Yi(/i; U), is in practice [/-independent and coincides with 
its 17 = limit X i(h; 0) ~ [A /h] 2 (such that X i{h; U) oc z(0)[A /h} 2 for fields h > f ). 

In summary, we have developed a local moment approach to static magnetic properties of 
the Anderson model. The approach encompasses all field and interaction strengths and, its 
undoubted simplicity notwithstanding, we find in particular for the strong coupling Kondo 
regime that it leads to good quantitative agreement with exact results obtained from sophisti- 
cated methods such as the Bethe ansatz P,[JT[. However the LMA can handle in addition both 
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non-intcgrable models and dynamical properties [6-8] . Its success in capturing static magnetic 
properties, particularly in the strong coupling regime, suggests its efficacy as a viable route 
to e.g. the field dependence of single-particle dynamics. We will discuss this important issue 
in a subsequent publication. 
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